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1 Studying the reason, why single-layer molybdenum disulfide (M0S2) appears to fall short of its 

promising potential in flexible nanoelectronics, we found that the nature of contacts plays a more 
f"^ ■ important role than the semiconductor itself. In order to understand the nature of MoS2/metal 

£SJ ' contacts, we performed ab initio density functional theory calculations for the geometry, bonding 

~. , and electronic structure of the contact region. We found that the most common contact metal (Au) 

is rather inefficient for electron injection into single-layer M0S2 and propose Ti as a representative 
example of suitable alternative electrode materials. 

Q\ ', PACS numbers: 73.40.Ns, 73.20.At, 73.22.-f, 81.05.Hd 

Contrary to popular perception, contacts often play a more crucial role in nanoscale electronics than the semicon- 
ducting material itself [3, ■ Whereas contacts in Si-based devices are no longer considered a problem after many 
decades of optimization, engineering optimum contacts to electronic nano-devices consisting of silicon^ or carbon 
(e.g. nanotubes or graphene) [J, 0] has become a major challenge for the field. More recently, the layered molyb- 
denum disulphide (M0S2) compound, which is structurally very flexible, has emerged as a promising alternative to 
silicon-based, carbon-based, and molecular electronics 0, Bulk M0S2, a well-established low-cost lubricant, has an 
indirect band gap of 1.2 eV0] and a rather high carrier mobi lity 3. In contrast to the bulk material, the observed 
electron mobility in single-layer M0S2 is unexpectedly low [HI. ITo|. 

Here we propose that the observed low electron mobility in M0S2 may not represent an intrinsic property of the 

■ semiconducting single layer, but was possibly biased by unfavorable contacts, which can dominate the electronic 
characteristics of MoS2-based nano-electronic devices. Our ab initio density functional theory (DFT) calculations for 
the electronic structure of MoS2/metal contacts indicate that Au, the most common contact metal in this system^], 
forms a tunnel barrier at the interface, which suppresses electron injection into M0S2. This is possibly the true reason, 
why the observed carrier mobility in single- layer M0S2 is lower than expected (El. flpj] . 

Searching for better contacts than provided by Au, we focusscd on metals with a low work function that would 
efficiently inject electrons into the conduction band of M0S2. Among transition metals with d orbitals that may 
favorably mix with the Mo4d states, we identified Sc, Ti and Zr as suitable candidates. Among these, Sc and Zr are 
. less suitable due to a large lattice mismatch, and Ti emerges as an ideal candidate with only 1% mismatch to M0S2. 
1 As we will show, the MoS2/Ti interface displays a much higher density of delocalized states at Ep than the M0S2/AU 
. contact. Similar to Au, Ti fulfills also other criteria required of a good contact material in electronics, such as high 
1 conductivity and chemical, thermal and electrical stability. Therefore, Ti is being used widely as a contact metal in 

■ modern microelectronics 
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Our DFT calculations use the Perdew-Burke-Ernzerhof form of the exchange-correlation functional as imple- 
mented in the SIESTA codeJT^]. A similar approach had been successfully used to characterize transition metal chalco- 
genide nanowires ID, 14 1 and their contacts to metal electrodes [lEj]. The behavior of valence electrons was described by 



norm-conserving Troullier-Martins pseudopotentials [1 61 ] with partial core corrections. We used a double-zeta basis, 
including initially unoccupied Mo5p orbitals. The Brillouin zone of the periodic array of MoS2/metal slabs, separated 
by a 33 A vacuum region, was sampled by a 8x16x1 k— point grid. We limited the range of the localized orbitals 
in such a way that the energy shift caused by their spatial confinement was no more than 140 meVfT^. The charge 
density and potentials were determined on a real-space grid with a mesh cutoff energy of 200 Ry, which was sufficient 
to achieve a total energy convergence of better than 0.1 meV/unit cell during the self-consistency iterations. 

The supercell geometry of the relaxed commensurate contact region between M0S2 and the close-packed surfaces of 
Au and Ti is shown in Figs.QJa) andQJb). We have represented the metal surfaces by 6-layer slabs. The single-layer 
M0S2 consists of a molybdenum monolayer sandwiched between two sulfur monolayers. In each unit cell we distinguish 
between two types of S atoms: two S atoms are located in the on-top site and six S atoms in the hollow site of the 
metal substrate. When optimizing this structure using the conjugate gradient technique, all atoms except the bottom 
four layers of the metal slabs were allowed to relax. In principle, the precise geometry at the metal-semiconductor 
interface may be affected by our choice of the exchange-correlation functional, which does not describe dispersive 
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FIG. 1: (Color online) Side view of the relaxed contact region at the interface between M0S2 and the (a) Au(lll) and (b) 
Ti(0001) surface, (c) Binding energy E per interface metal atom as a function of the separation d between between M0S2 and 
the Ti(0001) and Au(lll) surface. 



interactions accurately. Since the interface bonds are either semi-covalent or covalent and contain in our estimate 
only <20% van der Waals character, we expect the optimized geometry to be adequate for our study. We noticed 
that the relaxation within the M0S2 structure after contacting the metals was very small. 

The following major factors determine the electronic transparency of contacts: favorable interface geometry and 
bonding, the electronic density of states, and the potential barrier at the interface. Strong interconnects are especially 
important when contacting flexible semiconductors such as M0S2, which is known to form both planar and tubular 
nanostructures [3 ■ Favorable geometry precludes a small lattice mismatch at the interface, and should maximize 
the overlap between the states at both sides of the interface. The density of states (DOS) at the Fermi level (Ep) 
should be large throughout the interface region, forming delocalized states with low effective electron mass in order to 
efficiently transfer electrons between the metal and the semiconductor. The potential barrier at the interface should 
be as narrow and low as possible to maximize current injection. In the following, we analyze each of these factors for 
the interface between M0S2 and Au and Ti as contact metals. As will become clear in the following, Ti turns out to 
be superior to the commonly used Au as a contact metal. 

Contrary to chemically not saturated sulfur, which forms favorable thio bonds to Au, the sulfur in M0S2 is fully 
saturated and does not bond strongly to Au. This is reflected in the fact that the shortest distance of 2.62 A between 
S in the on-top site and the Au atoms directly underneath is about 0.2 A longer than the sum of the S and Au 
covalent radii. The distance between the S atoms in the hollow sites and their closest Au neighbors is significantly 
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FIG. 2: (Color online) Partial electronic density of Mo and S states, which are relevant to bonding and charge injection, in (a) 
the single-layer M0S2, (b) the M0S2/AU system, and (c) the MoS2/Ti system. Only a narrow energy region around the Fermi 
level Ef is shown. 

larger, namely 3.15 A. The average separation between the top layer of Au and the Mo layer is 4.21 A, large enough 
to suppress any efficient wavef unction overlap. 

When compared to Au, the equilibrium separation between Ti and M0S2 is much lower. In this case, the majority 
S atoms, which are in the hollow site, play the key role for adjusting the interlayer separation, trying to replicate the 
environment of Ti in the stable compound TiS2 . With the equilibrium distance between S atoms at the hollow site 
and its closest Ti neighbors of 2.54 A, the optimum separation between M0S2 and Ti is about 2.0 A, much shorter 
than the sum of the Ti and S covalent radii of 2.38 A. The resulting repulsion between the minority S atoms, which 
are in the on-top site, pushes away the Ti atoms underneath, reaching an equilibrium distance of 2.34 A, close to the 
sum of the respective covalent radii. The average Mo-Ti distance is 3.57 A, which is 0.64 A shorter than the Mo-Au 
distance, indicating favorable conditions for a large wavefunction overlap. 

To characterize the contact bond strength, we define the binding energy E between the metal and the M0S2 layer 
as the total energy difference between the combined and the isolated systems and display our results in Fig. QJc) . We 
find that the binding of M0S2 to Au is considerably weaker than to Ti, with the binding energy per surface metal 
atom of 0.36 eV in case of Au as compared to 0.98 eV in case of Ti. 

The electronic transparency of a contact can be quantified in a quantum transport calculation. For low bias voltages, 
a suitable approach may involve calculating the equilibrium Green's function, which to a large extent reflects the 
electronic density of states near Ep and the degree of derealization of these states within the contact region. For a 
detailed insight into the reason, why some contacts are better than others, we proceed with a careful analysis of these 
quantities. 

The DOS projection onto selected Mo and S orbitals is presented in Fig. [2] for the single-layer M0S2, M0S2/AU and 
M0S2/TL The bottom of the conduction band and the top of the valence band of the single- layer M0S2 is dominated 
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FIG. 3: (Color online) Electronic structure at the interface between M0S2 and (a-c) Au and (d-f) Ti. Contour plots of the 
charge density pi associated with states in the energy range Ef — 0.1 eV< E < Ef + 0.1 cV in planes normal to the interface 
in (a) M0S2/AU and (d) M0S2/TL Average value of < pi > (z) in planes parallel to the interface of (b) M0S2/AU and (e) 
M0S2/TL Average electrostatic potential <Vi> (2) in planes normal to the interface in (c) M0S2/AU and (f) M0S2/TL The 
dotted line in the panels indicates the location of the sulfur layer closest to the metal and the dashed line the position of the 
topmost metal layer. 



by Mo4g? z 2 states, with the other Mo states playing a minor role. Since the work function $(MoS2) = 5.2 cV is larger 
than that of most metals, the electronic transparency of the contact is maximized when electronic states at the Fermi 
level of the contact metal align and strongly overlap with the Mo4g? z 2 states near the bottom of the conduction band. 

As seen in Fig. [3Jb), upon making contact with Au, the MoS2/metal interface becomes metallic. The Fermi level 
of the combined system shifts upwards, to about 0.1 eV above the bottom of the conduction band of M0S2. The 
states near Ep display dominant Mo4d z 2 character with only a small admixture of S3sp states. Consequently, electron 
injection into the semiconductor will involve primarily the Mo4d z 2 states that, according to Fig. Hfb), display a low 
partial DOS at Ep. The corresponding low carrier density near Ep of the interface is depicted in Fig. [3ta) and 
Fig. [Ib). 

As seen by comparing Fig.[2Jb) for M0S2/AU and Fig.^c) for MoS2/Ti, Ti as contact metal modifies the electronic 
states near Ep much more than Au. In the MoS2/Ti system, the Fermi level is shifted upwards, to 0.25 eV above the 
bottom of the M0S2 conduction band. This is much higher than in the Au system and causes an increase in the DOS 
at Ep. The most striking difference to Au is a significant contribution of S3sp and Mo4d xy states near the Fermi 
level, which is associated with a strong S-Ti mixing. The contribution of Mo4.rf.cy states at Ep is nearly equal to that 
of the Mo4e? z 2 states. The broadening of the peaks in the DOS near Ep reflects an increase in the dispersion of the 
corresponding bands and suggests the formation of delocalized states at the interface. 

The character of the states, which determine the low-bias transport, is represented by the density p;(r) of the 
corresponding carriers in the left panels of Fig. (3{a) for M0S2/AU and Fig. [3Jd) for M0S2/TL As seen in Fig. EJa), 
which shows a detailed contour plot of pi, and Fig. [3jb), representing the average < pi > (z) in planes parallel to 
the interface, the charge carrier density in the interface region between M0S2 and Au is very low. Consequently, 
the electron transport across the M0S2/AU contact is mainly of tunneling nature. Since according to Fig. EJb) the 
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electron injection into the M0S2 layer proceeds exclusively via the Mo4<i z 2 states, the tunnel barrier from Au to M0S2 
is very wide. 

In striking contrast to the M0S2 / Au interface, the charge carrier density in the interface region between M0S2 and 
Ti is much larger. This is seen especially when comparing the averaged carrier density < pi > in planes parallel to 
the interface in Fig. ELb) for M0S2/AU and Fig. [21(e) for M0S2/TL Of particular interest is the difference between 
the electron density at the interfacial sulfur layer, denoted by the dotted line in Fig. [3] The nearly vanishing value 
of < pi > at this location in M0S2/AU increases by an order of magnitude in MoS2/Ti, thus turning tunneling into 
resonant transport. The carrier density derealization in MoS2/Ti, anticipated above due to the DOS broadening, 
corresponds to a metallization of the interface. This in turn enables direct charge injection into the M0S2 layer, for 
which the actual distance between Mo and Ti atoms becomes irrelevant. 

In order to complete the analysis of the contacts, we investigate the electrostatic potential regarding the existence 
of barriers at the metal-semiconductor interface and show the results in Fig. [3jc) for M0S2/AU and Fig. |3jf) for 
M0S2/TL Since we observe not only a net charge transfer across the metal-semiconductor interface, but also changes 
in the electronic structure due to the covalent interaction, these barriers are not ideal Schottky barriers, but rather 
more general contact tunnel barriers. We define the height of the contact tunnel barrier as the difference between the 
averaged potential at the top metal layer, indicated by the dashed line in Fig. [3J and the maximum of the averaged 
potential between the metal and the neighboring sulfur layer. 

As anticipated, the tunnel barrier at the interface between the semiconducting M0S2 layer and the Au surface, 
which is shown in Fig. Etc), is relatively high (1.03 eV) and wide. Since the S states of the bottom M0S2 layer are 
nearly absent near Ep in the M0S2/AU system, the true barrier is even wider than the 1.59 A value at half- height 
shown in Fig. Etc). In this case, tunneling involves direct charge transfer from Au to Mo states across two barriers. 

Electron injection from Ti to M0S2 is completely different. As seen in Fig. Etf), electrons in this system have 
to bypass the much lower (0.45 cV) and narrower (0.9 A) barrier to reach the delocalized states at the MoS2/Ti 
interface. In comparison to Au as contact metal, the significant reduction of the barriers at the interface with Ti will 
significantly improve the electronic transparency of the contact. 

Even though Au is commonly believed to be the ideal contact metal to many sulfur-terminated systems including 
multi- layer M0S2, our study shows that the opposite is true when contacting single-layer M0S2. A multi-layer system 
is preferentially contacted from the side, where Au can bond chemically to not saturated sulfur atoms at the edge. 
Since contacting a single-layer from the side is insufficient for good electron injection, the preferred geometry is a 
top contact discussed here. In this scenario, we identified unexpected qualitative differences between different contact 
metals in the way they inject carriers into M0S2. 

The basic difference is that between an inefficient tunnel contact in M0S2/AU and a low- resistance ohmic contact 
providing a direct injection channel in M0S2/TL We discuss MoS2/Ti only as a representative example of an optimum 
designer contact that is superior to the state-of-the-art. Good contact metal candidates must, of course, first fulfill 
macroscopic criteria such as high conductivity and chemical, thermal and electrical stability. Additional criteria for 
an optimum contact in nanoclcctronics, which we find fulfilled in the case of Ti, include a favorable interface geometry 
and bonding. In terms of electronic structure, an optimum contact has a high density of delocalized states across 
the interface at the Fermi level of the combined system, corresponding to a minimized or non-existent tunnel barrier 
between the two materials. 

In conclusion, we performed ab initio density functional theory calculations of M0S2/AU and MoS2/Ti contacts to 
study the reason, why single-layer molybdenum disulfide appears to fall short of its promising potential in flexible 
electronics according to recent experiments (EL fTo|] . We found that the nature of contacts plays a more important role 
in these systems than the semiconductor itself. Our calculations for the geometry, bonding and electronic structure 
of the contact region suggest that the most common contact metal (Au) forms a tunnel contact to single-layer M0S2 
and thus is rather inefficient for electron injection. We find that Ti is a suitable alternative as electrode material, 
since it forms a low-resistance ohmic contact. We also provide specific criteria for selecting materials besides Ti that 
should optimize the electronic transparency of the contact. Higher contact transparency reduces the required bias 
voltages for operation and may also improve the frequency response of these structurally flexible electronic devices, 
which may eventually open new horizons for electronics based on transition metal chalcogenides. 

DT was partly supported by the National Science Foundation Cooperative Agreement #EEC-0832785, titled 
"NSEC: Center for High-rate Nanomanufacturing" . GS was partly supported by the European Research Council 
(ERC - project INTIF 226639). IP was supported by the Science Foundation of Ireland (SFI) and CRANN. Com- 
putational resources for this project were provided by the ZIH Dresden and the Trinity College High Performance 
Computing Center (TCHPC). 



6 



* E-mail: tomanek@pa.msu.edu 
[1] F. Leonard and A. A. Talin, Phys. Rev. Lett. 97, 026804 (2006). 
[2] Y.-F. Lin and W.-B. Jian, Nano Lett. 8, 3146 (2008). 

[3] U. Landman, R. N. Barnett, A. G. Scherbakov, and P. Avouris, Phys. Rev. Lett. 85, 1958 (2000). 
[4] D. T. N. Nemec and G. Cuniberti, Phys. Rev. Lett. 96, 076802 (2006). 

[5] B. Radisavljevic, A. Radenovic, J. Brivio, V. Giacometti, and A. Kis, Nature Nanotech. 6, 147 (2011). 
[6] Y. Yoon, K. Ganapathi, and S. Salahuddin, Nano Lett. 11, 3768 (2011). 
[7] K. K. Kam and B. A. Parkinson, J. Phys. Chem. 86, 463 (1982). 
[8] R. Fivaz and E. Mooser, Phys. Rev. 163, 743 (1967). 

[9] V. Podzorov, M. E. Gershenson, C. Kloc, R. Zeis, and E. Bucher, Appl. Phys. Lett. 84, 3301 (2004). 
[10] K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. Khotkevich, S. V. Morozov, and A. K. Geim, Proc. Natl. Acad. 

Sci. USA 102, 10451 (2005). 
[11] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

[12] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Junquera, P. Ordejon, and D. Sanchez-Portal, J. Phys. Cond. Matter 
14, 2745 (2002). 

[13] I. Popov, T. Yang, S. Berber, G. Seifert, , and D. Tomanek, Phys. Rev. Lett. 99, 085503 (2007). 
[14] T. Yang, S. Okano, S. Berber, and D. Tomanek, Phys. Rev. Lett. 96, 125502 (2006). 

[15] I. Popov, A. Pecchia, S. Okano, N. Ranjan, A. D. Carlo, and G. Seifert, Appl. Phys. Lett. 93, 083115 (2008). 
[16] N. Troullier and J. L. Martins, Phys. Rev. B 43, 1993 (1991). 

[17] E. Artacho, D. Sanchez-Portal, P. Ordejon, A. Garcia, and J. M. Soler, Phys. Stat. Sol. (b) 215, 809 (1999). 
[18] R. Tenne, L. Margulis, M. Genut, and G. Hodes, Nature (London) 360, 444 (1992). 



